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CTTROBPCTIOff 

Optimisation of structural members has been a very intriguing 

r 

topic of investigation in the past two decades References 1-2).. 
Development of new structural, materials, such as composites, and a 
great need for light weight structures has made it even more important 
to find optimum designs using composite materials. 

There are several different levels of abstraction at which the 
basic structural design problem can be approached. The most common one 
is to consider the optimum design of structural elements such as columns 
and plates and composite structures such as box beans and panels for 
prescribed loads and prescribed overall (leading) dimensions. Thus for 
individual elements, the optimum design analyses result in the specifi- 
cation of the cross sectional dimensions for a given loading index. In 
certain applications, it is meaningful to relax one of the leading 
dimensions of a composite structure to find a design of absolute minimum 

weight. In this case Eecker has obtained results for the optimum 

t 

structural chord for a box bean (Reference 3) and an optimum diameter 

f, 

for a cylinder in bending (Reference U) by essentially adjusting the 
value of the loading index. 

The concept of loading indices and efficiency factors have been 
proved very useful for the conventional isotropic materials. The 

r 

development of these concepts is attributed to Zahorski (Reference 5) • 


1 
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They have been used very effectively by Farrar, Shanley, Gerard and 
f 

others (References 6, T» 8), The loading index concept is applied in a 
minimum veight or efficiency analysis by expressing the quantity to be 
minimized (weight) or maximized (stress) in terns of the prescribed 
dimensions and loads. In doing so, the general approach used by many 

ti 

investigators i-s to reduce the number of unknown dimensions to two or 
three by making suitable guess so as to the ratio of various dimensions 
in order to get a closed form solution. 

In the present study a relatively new approach of structural, 
optimization has been used to optimize the weight of a simply supported, 
corrugated hat stiffened composite panel under uniaxial compression. 

This approach consists of the employment of nonlinear mathematical 
programming techniques to reach an optimum solution. This approach is 
in contrast to the one for which a closed form is attempted, since for 

t 

the later no simplifying assumptions are required in general with regard 

* 

to the cross-sectional dimensions. However in the present work some 
simplifying assumptions in the stress analysis are made to effect faster 
convergence to an optimum solution. With these simplifying assumptions 
the number of unknown design parameters is reduced to twelve for the 
purpose of optimization. Since the number of unknown parameters in the 
present problem is twelve as compared to two or three in the loading 
index approach, twelve simultaneous equations are needed to get values 
of all the unknown parameters. Hence, in the loading index approach, 
either further simplifying assumptions have to be made for the dimen- 
sions of the cress section or a more involved stress analysis describing 
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the behavior of failure is required. For examples in buckling problems 
additional inodes of failure have to be considered. 

In the present analysis, a computer code (Reference 9) called 
AESOP (Automated Engineering and Scientific Optimization Program) is 
used for the optimization studies, 

AESOP consists of several optimum search algorithms. Depending 
upon the behavior of the performance function (weight), one or a 
combination of search algorithms can be used to find the parameters 
(design variables) which will minimize the performance function. 

AESOP is used to optimize the design parameters of the panel. Then, as 
a check on the effect of the simplifying assumption, the critical load 
for the optimized panel is obtained from EUCLASP-2 (Reference 10) and 
compared with the specified panel loads. Good correlation was obtained. 

Unfortunately, no optimization results are available for all 
composite panels, for the pur-pose of comparison. The results for alumi- 
num panels are available (Reference 2 ). A comparison of the present 
results was made with the available results and good correlation was 


found. 
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GENERAL APPROACH OF STRUCTURAL OPTIMIZATION 

Portions of this chapter closely follow the material contained in 
Reference (9). ■ 

2.1 Concept of optimization 

In general, any optimization problem can be thought of as minimi - 
zation or maximization of a performance function. For example, in 
structural problems, weight and stress are the performance functions. 
Similarly in rocket design, the range of the rocket may be taken to be a 
performance function. In all such optimization problems the ultimate 
aim is to find the value of design parameters which will optimize the 
performance function. 

2.2 Optimization approaches 

There have been several different approaches used by many investi- 
gators to reach to an optimum solution. One of the most powerful 
techniques used until recent years for the optimization of conventional 
structural members is the loading index approach (Reference 2). In 
this approach, the loading index is expressed in terms of weight or 
stress and the dimensions of the structural member. In general this 
can be written as 

loading index = efficiency factor x (weight index) 11 (l) 


k 
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la this equation the efficiency factor is a function of geometric 
properties of the structure under consideration. These may he treated 
as independent variables in minimization of the veight. The weight 
index is designated as the nondimens ional weight function. Minimization 
of weight is achieved by maximization of the efficiency factor. In 
equation (l) n is an exponent, whose value depends on the structure 
under consideration. 

The basic assumption used in arriving at equation (l) for a given 
structural element of buckling problem is that, for optimum design, 
at least two lowest modes of instability are simultaneously critical under 
the applied loading (References 7, 11). It should be noted that in this 
approach the number of unknown parameters is reduced to two or three 
in order to get a closed form solution. The reason for doing so is 
obvious. For more unknown parameters, the problem becomes more complex 
and it becomes impossible to get closed form solutions because more 
modes of failure have to be considered in order to get additional equa- 
tions to determine the unknown parameters. 

In the present problem the loading index approach can not be used 
effectively because the number of unknown parameters is large. A more 
recent approach of structural optimization is to use nonlinear mathe- 
matical programming techniques. In this approach we are concerned with 
the maximization or minimization of a pay-off or performance function <j> 
of the form 

4> ~ Q(<* ± ) » 1 ~ 1» 2, — — — , n (2) 
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subject to the array of constraints 

C j S W » J = 2, , P (3) 

The cu are the independent design variables whose values are to be 

determined so as to maximize or minimize the performance function 4>(ct^)» 

subject to the constraints of equation (3). The ou may be thought 

- N 

upon as the components of a control vector, a, in the space R of 
dimension E. Since maximization of a function is equivalent to mini- 
mization with a change of sign, it will be sufficient to discuss the 
case in which performance function is to be minimized. 

Multivariable optimization problems involving inequality regional 
constraints relating the design variables may also be "encountered as 
follows 


00 

The inequality constraints define a region of the control space within 

which the solution must lie. For example, in structural buckling 
problems, if the design variables are taken to be cross-sectional 
dimensions, then these dimensions neither can be less than or equal to 
zero nor become infinitely large. So the above limits bound the region 
in which these variables must lie. 

Inequality constraints on the functions of independent variables 
similarly restrict the region in which the optimal solution is to be 



? 


obtained. In this case 


£ E K (a i> 


( 5 ) 


For example, in structural buckling problem, various modes of buckling 
* 

modes will constitute constraints, such that the structure is capable 
of carrying the design load. These constraints will be function of 
independent design variables. 

Inequality constraints can be used to restrict the search region 

directly, or, alternatively, they may be transformed into equality 

constraints. Several transformations may be used for this purpose. 

For example, let an equality constraint, C„, be defined by the 

ft. 

transformation 


(p k- f / ■' f k < p £ 


°K 


P K i F K i F K 




<4- f k )2 • 4 < f k 


( 6 ) 


constraint C„ to zero will result in the constraint of equation 
ft 

(5) being satisfied. 

Problems involving equality constraints can be treated as un- 
constrained problem by replacing the actual performance function, 
<j>(a^), by a penalized performance function, where 
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P 

<j>* = 4> + Z U C (7) 

5=1 3 3 

It can "be shown that, provided the positive weighting multipliers 
are sufficiently large in magnitude, minimization of the performance 
function subject to the constraints of equation (3) is equivalent to 
minimization of the unconstrained penalized performance function defined 
by equation (7). This approach permits search techniques for finding 
unconstrained minima to he applied in the solution of constrained minima 
problem at the cost of some increased complexity in the behavior of the 
performance function. The weighting multipliers are determined 

adaptively on the basis of response surface behavior. 

Alternatives to this approach are available, notably Bryson’s 
approach to the steepest-descent search (Reference 12). This method 
has been exploited in connection with the numerical solution of varia- 
tional problems encountered in the optimization of aerospace vehicle 
flight paths (Reference 13). However the use of such techniques implies 
smoothness of the response surface. This smoothness can not be assumed 
in the problem of structural optimization in general; hence, the less 
restrictive penality function approach of equation (7) is used. A 
detailed discussion of solution techniques is presented in appendix A, 
however, some of the search algorithms will be discussed below. 
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2.2.1 Adaptive creep search 

This search is a form of small scale sectioning; however instead 
of locating the position of the one-dimensional extremal on each section 
parallel to a coordinate axis, the coordinate is merely perturbed by 
small amount, Au^, in the descending direction. 

The search commences with a small perturbation in one of the 
independent variables, ct^; a positive perturbation is first made; if 
this fails to produce a performance improvement, then a negative pertur- 
bation is tried. If neither of the perturbations produces an improved 
performance value, the variable retains its original value, and Aa^ is 
halved. If a favorable perturbation is found, the variable is set 

to this value , and Act^ is doubled. The process is repeated for each 
independent variable in turn, the order in which the variables are 
perturbed being chosen randomly. At this point an adaptive search 
cycle is complete, and the cycle is then repeated. A two-dimensional 
illustration of this search is presented in figure (l). In the 
particular problem illustrated, the method converges rapidly reaching 
the neighborhood of the extremal within six evaluations. 

The search algorithm can be written in the form 

(S„ - T ) 

Aa r = 2.0 r r . (DP) (8) 

where S is the number of cycles in which the search has successfully 
th 

perturbed the r independent variable, and is the number of 

*fc$l 

cycles in which the perturbation of the r 


variable has proved 
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unsuccessful. Here, the scalar quantity (DP) merely defines an initial 
perturbation for each independent variable. Once started the search 
proceeds inevitably to its conclusion, the perturbation in each inde- 
pendent variable being adaptively determined according to equation (8) 
on the basis of the performance function response contour behavior 
encountered during the particular problem solution. This search can be 
quite efficient vhen used in combination with the pattern search 
acceleration procedure. 

2.2.2 Pattern search 

In the present work, pattern search refers to a search vhich 
exploits a gross direction revealed by one of the other searches. The 
search algorithm is 

Act.. = (a* - c£) • (DP), i * 1, 2, ...... N (9) 

2 1 

where and are the components of the control vector before 

and after the use of a preceding search technique. This is illus- 
trated in figure (l) following an adaptive search. The combination of 
an adaptive search and a pattern search in the problem illustrated 
leads directly to the neighborhood of the extremal. Repeated adaptive 
search on the other hand, would be a very slowly converging process 
due to the orientation of the contours with respect to the axes of the 
independent variables. It may be noted that a simple rotation of the 
independent variable axes by ^5° results in adaptive creep alone 
becoming a rapidly converging process in this example. The present 
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discussion of optimization concept is rather superficial. Detailed 
treatments may he found in (References 9 S 22). 

2.3 General structural optimization cycle 

Figure (2) shows a typical optimization cycle. First of all the 
geometry (e.g. flat panel with corrugated hat stiffeners) of the struc- 
ture, the loads, and the material are specified. An attempt is then 
made to find the values of design variables, which will minimize the 
weight of the structure. Figure (2) shows some of the design variables 
for a composite panel subjected to compressive loads. These design 
variables are 

1. Cross-sectional dimensions : In the optimization problem one 

has to find the dimensions of the cross-section which will minimize 

* 

the weight of the panel. 

2. Filament orientation : This means, the orientation of filaments 

with respect to a reference axis. 

3. Percentage of different orientations : This simply means the 

percentage of differently oriented laminates required in a Lamina . 

For example, if only two kinds of laminate orientations are used, 
say 0° and 90° one has to know that, how much of each is needed in 
a structural member to obtain the least weight design. 

k. Laminating sequence ; This means the sequence in which differently 
oriented laminates are arranged to make a structural member. 

There could be more design variables depending upon the need of the 
problem under investigation. In order to find the value of these design 
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variables one has to iterate several times in such a way that these 
values represent the design of minimum weight. We note that a complete 
stress analysis of the structure has to be made during each iteration 
cycle. In order to cut down on the iteration time, it will be very 
necessary to make some simplifying assumption for the purpose of 
stress analysis. This is done in the next chapter in which a detailed 
formulation of the problem is given. 



Ill 

ANALYSIS 

3 * I Formulation of the problem 

The problem under consideration here, is that of optimization of 
a simply supported all composite corrugated hat stiffened panel under 
uniaxial compression. Figure (3) shows the panel under consideration. 
The material used in the analysis is graphite/epoxy. 

3.1.1 Basic assumptions 

I. All the panel members are thin plates simply supported on all 
four edges. 

?., All the panel members are orthotropic and have constant 
thickness. 

3. Only three kinds of laminate orientations are used, namely 0°, 
+45° * “45° , relative to the axial direction. 

4. Yield strain in compression for any panel member is equal to 

yield strain of 0° laminates irrespective of the percentage of 0° 
and laminates. 

5. Each panel member is assumed to have only three layers. 

6. Laminate layup in each panel member is assumed to be +45°, 

0°, +45°. 

J. Effect of Poisson’s ratio is neglected in calculating the 
load carried by each panel member. 
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8. Panel is assumed to behave like a wide column for the purpose 
of Euler buckling analysis, 

9 - Torsional and local cripping failures of the panel modes are 
ignored. 

3.1.2 Performance function 

Under the assumption of vide column behavior only one pitch of 
the stiffner spacing is required for the purpose of further analysis . 
Figure (4) shows a representative cross-section of the panel. 

The performance function in cur analysis is a function of weight. 

It is chosen to be the weight per unit area per unit width of the panel. 

$ = ££ do) 


3.1.3 Design Variables 

Taking into consideration assumption (3) » (5)> and (6), the number 
of unknown design variables is reduced to twelve. They are (also see 
Figure 4) 


1. Width of each panel member, b^ 

2. Thickness of each panel member, t^ 

3. Percentage of +45° laminates in each panel member, f^ 
where i = 1,2, 3,4 


(H) 


3,1. A Constraints : 

The panel must meet certain failure criterian and practical 



15 


requirements in order to be a valid design. They are as follows 

1. Local buckling load of each panel member should be greater 
than or equal to the applied load. 

P, > P„ (12) 

4 1 a l 

2. Euler buckling load of the total panel should be greater 
than or equal to the applied loading, 

N > N (13) 

jSuler a 

3. Applied strain of the total panel should be less than or 
equal to cutoff, or yield strain. 

e < e (lh ) 

h. Stiffner spacing should be greater than or equal to b^ 

(Figure )j) 

b>b 3 (15) 

5« Value of design variables should be limited in a region of 
practical interest. 



16 


»i<\< b= 


t£ < t t < tj 


f L < f < f® 




J 


(16) 


For example, percentage of +U5° laminates (f^) can not be less than 
zero and greater than hundred. 

The value of all the above mentioned parameters are obtained through 
the use of a simplified stress analysis, as discussed below. 


3.2 Stress Analysis 

3.2.1 Load in each nanel member 

Lets assume N is the load intensity per unit width, o . is the 

X &1 

axial stress in each panel member, and P is the total load per unit 
stiffner spacing. Nov we can write 


and 


P = N • b 
x 


(17) 


P = 


4 

t 

i=l 



(18) 


Because of compatibility considerations, the strain in each panel 
member has to be equal. Here we will neglect the effect of Poisson’s 


ratio . Hence 
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o , 

e * —■ i = 1, 2, 3, U (19) 

E i 

From Equations (18) and (19 ) ve get 

k 

P = £ e E. A. 

1*1 1 1 

Solving this equation for the strain e and using Equation (17) » we 
obtain 


e 


N b 
x 


£ E. A. 
i=l 


( 20 ) 


Finally the load P in each panel member is then given by 

a i 



N x b F i A i 
~k 

. E F i A i 
x=l 


( 21 ) 


3.2.2 'Local buckling 

Each panel member is assumed to be orthotropic and simply supported 
on all four edges. Hence from orthotropic plate theory (reference 13 )> 
we can write 


and 


- fr (i ^T^ + + b 66> 

D % 



? r^ bt 


(22a) 


(22b) 


where , 
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cr^ - is the local buckling stress 
b - is the width of the plate 
t - is the thickness of the plate 

- are the bending stiffness coefficients 
“ is the local buckling load 

Here for the sake of simplicity subscript i has been, omitted, 
this equation applies to each panel member. 


However 


3.2.3 Euler buckling 

e 

We will consider one pitch of stiffner spacing for the purpose of 

4 

Euler buckling analysis. Since each panel member can have different 
percentage of + 1*5° laminates, each panel member is liable to have 


different value of Young's modulus. In order to find the Euler buckling 


load, we will use the equivalent area approach to find the effective 
Young's modulus and effective area of each panel member. Then we can 
employ 


JT 


X Eul' 


er 




b L 


(23) 


In this expression E^ * I is the effective stiffness. Next we obtain I. 
Let us asstime 

y^ - distance of center of gravity of i^* 1 panel member from reference 
axis (shown in Figure It) 

Ej - Young's modulus for the 1 panel member 



19 


A^ - Area of the panel member 

y “ distance of the effective center of gravity of EA 
distribution from reference axis 
so ve can write 


X E A y 
i=l 1 1 X 


Z E. A. 
i=l 1 


By dividing both the denominator and the numerator by E^, we obtain 


* h h 

E ir iy i 


y . Ifl ^ 
y 4 E. A. 

s "t — 

i=i 


Let 


E^ 

3= '■ - 

i E, 


(24) 


therefore 


y = 


Z A* E 
1=1 


Z A* 
1=1 1 


From Figure (4) we see that 


(25) 


+ + 2 v + v Z 1/2 

b 2 = {{b - ( -A— ■ * )) + ( ■- 1 - 2 -- 3 )} 3 


(26) 


and 
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Hence 


" Vi 

A* = 2 V 2 V E I 
A| - 

A S - 2 Vh + V V E i 


(27) 


- *2 b 2/2 * A * b 2 + *£<*1 * b l t )/ 2 

y “ (A» + A* + A* + Ajj) 


( 28 ) 


In order to find the moment of inertia about the center of gravity 
of the cross-section, ve should adjust the effective vidth and the effective 
thickness of each panel member in such a vay that the respective distance 
of the center of the gravity of each panel member should remain unaffected. 
For example, with reference to the figure (5) we note the width of the 
horizontal panel member has been modified, but not the thickness. 

However, for the inclined panel members, the thickness will be modified. 
Hence ve can define 


b 3 ■ A 3 /b 3 


b 5 * A &/H 


(29) 


t* = a»/E 2 
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Therefore, the moment of inertia I, is given by 

1 “ Vl/12 + b ! * 4/6 + *3 * 3/12 

" *2 (b 2 - (t l + t l t >/2 }2/12 + A 2 (7 - b 2/2 )2 
+ A*(y - b 2 ) 2 + A* y 2 + Ajj{y - (t 1 + tg)/2} 2 (30) 

3.3 Discussion of minimization procedure 

At this point it is noted that the standard weight strength parameters 
are N^/L and the weight per unit area per unit length (W/bL ) (Reference 
7). The length of the panel is not by itself a design parameter. In 
the present analysis, we will assume that the loading 'N and length L, 
of the panel are known design parameters. However it will be shown in the 
optimization process that this approach will lead to the same weight 
strength plot. That is, N /!» is indeed the pertinent parameter, and not 
H x and L separately. 

In order to minimize the weight function of equation (10) , it will be 
necessary that all the constraints of equation (12) through (l6), be 
satisfied. After satisfying all the constraints and reaching a minimum 
solution, the outcome of the analysis will be the extremiying values of 
design variables of equation (11 ). For the purpose of reaching a minimum 
solution, an optimization computer program AESOP is used (see appendix 
A and Reference 9). 

It is very important to check the effects of the simplifying 
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assumptions. For this purpose another existing computer program BUCLASP- 
2 is used. (See appendix B and Reference 10). This program is capable 
of performing the buckling analysis of a biaxially loaded composite panel. 
A short description of the assumptions made in BUCLASP-2 analysis and the 
mathematical model required for the purpose of analysis is presented in 
Appendix-B. It should be noted at this point that BUCLASP-2 is not used 
as an optimization program, but it is used to predict the buckling loads 
of aD optimized panel. 
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NUMERICAL RESULTS 

In the present work following two cases of hat stiffened panels 
under uniaxial compression have "been optimized using "AESOP" . 

1. All-aluminum panel 

2. All-composite (Graphite/Epoxy) panel 

The material properties used for the purpose of analysis are presented 
in Table (l). 

•Tables (2) through (7) give the values of optimized design variables 
for various loading conditions and length of minimum weight Graphite/Epoxy 
panel. The values of the design variables for various loadings and length 
of the minimum weight aluminum panels are shown in Table (8). 

Tables (l) through (3) show the values of weight per unit area per 
unit length (W/bL ) and axial load per unit length per unit width (N /L) 

X 

for different lengths of the panel. It can be seen that for same value 
of N^/L the corresponding value of W/bL is similar in all three cases. 
This proves that and L do not have to be considered separately 

but only N^/L should be considered while obtaining these plots. 

Figure (6) shows a standard weight strength plot for a composite 
and an all aluminum panel under uniaxial compression. In Figure (6), 

N^/L is the load per unit width per unit length, and W/bL represents 
the weight of the panel per unit area per unit length. 
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Results of Reference (2) and the results obtained through the use of 

\ 

BUCLASP-2 (Reference 10) are also presented in Figure (6) for the purpose 
of comparison* For the all aluminum panel, a comparison was made with the 
results obtained by Crawford and Burns (Reference 2). As may be seen from 
figure (6), the present results show a slight weight advantage over the 
results obtained in Reference (2). This may be due to the fact that in the 
present analysis no assumption “was made with respect to the cross- 
sectional dimensions of the panel. It is also noted that Reference (2) 
uses the condition that for minimum weight, the local buckling stress in 
each panel member is set equal to the Euler buckling stress of the whole 
panel. In the present analysis no such condition for minimum weight design 

vas imposed. However it is interesting to note that the results obtained 

« 

after minimization process shoved that in fact for minimum weight, local 
buckling and Euler buckling stresses should be equal in each 
panel member. 

For all composite panels an attempt was made to make a comparison 
with available optimization results. Unfortunately, the author was unable 
to find such results. So instead, a comparison was made to study the 
effects of simplifying assumption employed in the present stress analysis. 
This was dene by determining the buckling load for the optimum panel by 
using the BUCLASP-2 computer program, which is devoid of such assumptions . 
A comparison vas then made between the buckling load obtained by BUCLASP-2 
for the optimum panel and the specified load, that was used in the present 
analysis to obtain the optimum panel. Figure (6) shows good correlation 
between the present results and these obtained through the use of BUCLASP-2. 
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The advantage of employing simplified stress analysis is, that it results 
in very s m a l l computational time. For example, vith the use of the 
simplified stress analysis, the runtime for 1500 iteration is about four 
seconds on the CDC-6600 computer. It is interesting to note that even 
for all composite panels, the results obtained through the use of BUCLASP-2 
show that the local buckling load and Euler buckling load of optimized 
panels is very close to each other. 

Examination of Figure (6) reveals that all composite panels weigh 
approximately half as much as all aluminum panels. This result for all 
composite panels is very useful for modern aircraft technology. It is 
hoped that the results of this study will lead to further investigation 
in the use of composite materials for various design problems. 

Figure (7a) and (7b) show that there are two different design possible 
for the same loading condition Figure (7a) pertains to a lightly loaded 
panel (II^/L ® 50) whereas figure (7b) applies to a heavily loaded panel 
(N^/L - 500). For each loading case, both of the designs weigh almost the 
same (see Tables 2 and 5)» but both have different values of design 
variables. This phenomenon allows for more flexibility during the design 
process and less weight penalities will be felt if practical constraints 
(e.g. manufacturing restrictions) are imposed on such panels. However 
one should make sure in such cases of multiple optimum designs, that 
these designs are not the result of the various assumptions made during 
the stress analysis. Figure (8) shows the buckling mode shapes for the 
two highly loaded panels (Figure 7b) . These mode shapes were obtained 
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from BUCLASP-2 , The panel in Figure (8a) is very deep as compared to its 
width and fails in a torsional mode. Since the torsional mode of failure 
was neglected in the stress analysis of the given panel, this, .panel is not 
a valid design, and can he ignored. The panel in Figure (8b) fails in a 
local buckling inode, such that all the panel members behave almost simply 

i 

supported, which is in accord with one of the simplifying assumption. 

We note the panel on the bottom of Figure (T) have one hundred 
percentage +1*5° Laminates in the skin and in the inclined panel members. 

f 

Also the thickness of the inclined members and skin is very small compared 

to the thickness of other panel members (See Tables 2-7). This suggests 

that most of the load is carried by 0° filaments, which is desirable in 

order to have most efficient panel. At this point it should be noted 

* • 

that under such condition neglection of the' effect of Poisson's ratio is 
a good assumption. Recall that this assvuaptidh was made in the calcu- 
lation of the axial load carried by each panel member. These panels 
also verify a very useful concept of reinforcing hat stiffened metallic 
panels. In reinforcing such panels strong load carrying material is added 
in the direction of the loading. For example, the reinforcement is added 
along the flanges and the skin connections in case of hat stiffened 
panels. 

Finally, it is noted that some of the assumptions made in the analysis 
of these panels did not effect the results to any significant amount. 

These assumptions are discussed next. 

1. The yield strain assumed to be that of all 0° filaments is 
considered to be a good assumption because most of the load is carried by 
0° filaments. 
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2. Laminating layup in each panel member turned out to be of no 
importance, because all the panel members have only one kind of laminate 
orientation, i.e, +4 5° or 0°, and never have both 0° and +45° filaments 
orientation. 
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CONCLUDING RELARKS 

In the present vork an attempt has "been made to discover some of the 
nev concepts in the optimum design of structural members, composed of 
composite materials . Since there are no available results for the purpose 
of comparison, it will be very desirable to carry out experimental verifi- 
cation of the results obtained in the present a n alysis. It should be 
noted that the results obtained in the present analysis are optimums but 
are not very practical. For example, it is not desirable to have all 0° 
filaments in any of the panel members. The imposition of such practical 
constraints will result in a heavier panel. If instead, however, the design 
variables are chosen in such a way that an •optimum design includes the 
practical constraints, relatively lighter practical designs may be found. 

For the case considered herein, the present analysis shows that 
composite panels are approximately twice as light as all aluminum panels . 

It is hoped that this result will inspire further investigations into 
the use of composites for optimum designs. In the present analysis some 
of the assumptions were very crude and need to be modified. For example, 
the yield strain for the whole panel was assumed to be equal to yield 
strain of all 0° filaments. This assumption in the present analysis did 
not effect the design to any significant amount, because most of the load 
is carried by 0° filaments. However this will not be true in general. 

This suggests a more general yield criteria is required. Furthermore it 
will be interesting to investigate panels with different geometry, boundary 
conditions and loading conditions. 
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APPENDIX A 


Optimization Technique 

This appendix is presented here for the sake of completeness. The 
contents of this appendix are available in reference (9). In this appendix 
techniques available in AESOP, for the solution of non-linear multivariable 
optimization problems are discussed. A wide variety of search algorithms 
have been devised for the solution of multivariable optimization problems. 
Many of these algorithms are restricted to the solution of linear or 
quadratic problems. Algorithms of this type must be supplemented by more 
general search procedures if generality of solution is sought. This is 
because engineering problems tend to lead to non-linear formulation with 
the possibility of discontinuities in both the performance function response 
surface and its derivative. Most of the searches which prove effective in 
these problems combine a direction generating algorithm, such as steepest- 
descent, with a one-dimensional search. Distance traversed through the 
control space in the selected direction is measured by a step-size, or 
perturbation parameter DP. The object of the one-dimensional search 
is to determine the value of DP which minimizes the performance function 
along the chosen ray and to establish the corresponding control vector. 

In practice, the diverse nature of non-linear multi-variable opti- 
mization problems leads to the conclusion that no one search algorithm 
can be uniquilly described as being the ’’best" in all the situations 
which may be encountered. Rather, a combination of searches, some of 
which may be of quite elementary nature, provides the most reliable and 
economical convergence to the optimal solution. 
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One dimens zonal search : Multivariable search problems are reduced 

to one- dimensional problems whenever a search algorithm is used to 
establish a one-to-one correspondence between the control vector and a 
single scalar perturbation parameter, DP. In such a situation 

a ± = (DP ) , i = 1, 2, N (Al) 

so that equation (2) becomes 

4> - 4>(a ± ) = 4>(DP) (A2) 

Similarly, the right hand sides of equation (3) and (7) become functions 
of the scalar perturbation parameter. 

The relationship, equation (Al), specifies a ray through the control 
space. As noted above, the objective of the one-dimer.sicnal search along 
this ray is to locate the value of DP which provides the minimum perfor- 
mance function value. 

Kumerie&l search for the one-dimensional minima can be carried out 
in a local fashion, by the l.'evton-Raphson method, for example, or by a 
global, search of the ray throughout the feasible region. The localized 
polynomial approximation is appropriate to the terminal convergence phase 
in a problem solution when some knowledge of the extremal 1 s position has 
been accumulated by the preceding portion of the search and the problem 
involves a smooth function. The global search can be used to advantage 
in the opening moves of a search. In the early phase of a search the 
object is to isolate the approximate neighborhood of the minimum performance 
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function value as rapidly as possible, usually with little or no fore- 
knowledge of the performance function behavior. One measure of the 
effectiveness of a search algorithm in such a situation is the number of 
evaluations required to locate the minimum point to some prespecified 
accuracy.' It can be shown that the most effective method of locating the 
minimum point of a general unimodal function is a Fibonacci search 
(reference 21). In this method, the accuracy to which the minimum is 
to be located along the perturbation parameter axis must be selected 
prior to the commencement of the search. Since the accuracy required is 
highly dependent of the behavior of the performance function, this quan- 
tity is difficult to prespecify. 

Prespeeification of the accuracy to which the extremal's position is 
to be located can be avoided for little loss in search efficiency by use 
of an alternative search based on the so-called golden section. (reference 
21). This is the method employed in the AESOP code one-dimensional search 
procedure. Search by the golden section commences with the evaluation of 
the performance function at each end of the search interval and at 
G - 2/(1 + /5) of the search interval from both of these bounding points. 
This Is illustrated in figure (Al) . 

The boundary point furthest from the lowest resulting performance 
function value is discarded. The three remaining points are retained, 
and the search continues in a region which is diminished in size by G. 

The internal point at which the performance function is known in the 
reduced interval will be at a distance G of the reduced interval from the 
remaining bounding point of the original interval for (l-G) = G . The 
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search can, therefore, be continued in the reduced interval vith a single 
additional evaluation of the performance function. It follows after Q 
evaluations of the performance function that the position of the extremal 
point will be known within R of the original search region where 


R = 


(A3) 


To reduce the interval of uncertainty to .00001 of the original 
search interval, about 27 evaluations of the performance function are 
required. For a reasonable number of evaluations of the performance 
function, this type of search is almost as efficient as a Fibonacci 
search. 

It should be noted that search by the golden section proceeds under 
the assumption of uniraodality; hence it will often fail to detect the 
presence of more than one minimum when the performance function is multi- 
modal. If more than one minimum does exist, the one located depends on 
performance behavior within the original search interval. 

Multiple Extremals on a One-Dimens icnal Ray : The one-dimensional 

section search described above is unable to distinguish one local 
extremal, from another; it will merely find one local extremal. This 
difficulty can be largely eliminated by the addition of some logic to 
the search, at least for moderately well behaved performance functions; 
that is, for functions having a limited number of extremals in the 
control space region of interest. An effective method for detecting 
multiple extremals is to combine the one-dimensional search with a 
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random one-dimensional serach on the same ray through the control space . 
This is illustrated in figures A 2 and A3. In Figure A2 the response 
contours of a performance function having two minima are illustrated 
together with the initial points used in a global one-dimensional search 
by the golden section method. The behavior of the function at these 
points is shown in figure Ak A3. The left hand minimum is not apparent 
from these points. If a single random point is added in the interval 
Lq» the probability of this point revealing the presence of the second 
minimum is 


P 1 = L X /L Q (Ah) 

for any point in the interval AB indicates the presence of a local 
minimum somewhere in the interval AB, and any point in the interval BC 
indicates the presence of a local maximum somewhere in the internal BC. 
In this latter case, there must be a minimum of the function both to the 
left and to the right of the newly introduced point. 

If random uniformly distributed points are added in the interval L^, 
the probability of locating the presence of the second minimum becomes 

P R = 1.0 - (1.0 - l x /l 0 ) r (A5) 

The function (L^yL^) is a measure of the performance function 
behavior. For a given value of this behavior function the number of 
random points which must be added to the one-dinensional search to provide 
a given probability of locating a second minimum can be determined. 
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The presence of multiple minima on a one-dimensional cut through an 
H-dinensional space does not necessarily indicate that the performance 
function possesses more than one minimum in a multi-dimensional sense. 

It may "be that the performance function is merely non-convex. This is 
illustrated "by figure Ah. The performance function behavior on the one- 
dimensional search in figures A2 and Ah is identical. In figure A2 this 
indicates the presence of tvo local extremals; in figure AU, a non-convex 
performance function. 

When a one-dimens icnal search detects the presence of multiple 
extremals in the local sense above, a decision must be made as to which 
of the apparent extremals is to be pursued during the remainder of the 
search. Here, without foreknowledge of the performance function behavior, 
logic must suffice. Typically , the left or right hand extremal, the 
extremal which results in the best performance, or even a random choice 
may be made. 

It should be noted that logic of this type is not currently 
available in the AESOP code. The AESOP one-dimensional search procedure 
has three distinctive phases. First, each search algorithm defines an 
initial perturbation using either past perturbation stepsi 2 e informa- 
tion or a perturbation magnitude prediction as in the quadratic search 
(Reference 9)* Second, a perturbation stepsize doubling procedure is 
employed until a point exhibiting diminishing performance is generated. 
Third, having coarsely defined the one-dimensional extremal position from 
steps one and/or two, a golden section search is employed to locate the 
extremal with reasonable precision. 
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Multiple extremals - general procedure ; The multiple extremal search 
technique included in AESOP is based on topologically invariant warping 
of the performance response surface. The response surface is warped in 
a manner which retains all the surface extremals hut alters their relative 
locations and regions of influence. The regions of influence of an 
extremal is defined as the hull or collection of all points which lead to 
the extremal if a gradient path is followed. Reducing the region of 
influence of an extremal decreases the probability of locating a point in 
the neighborhood of the extremal if points are chosen at random. Again, 
in an organized multivariable search, the probability of locating an 
extremal having a small region of influence is less than that of locating 
an extremal having a large region of influence. For example, suppose the 
extremals of the one-dimensional function of figure A5 -are to he determined 
in the range < a by the sectioning approach. The four initial 
values employed in this technique are denoted by f^ to f^. 

Following evaluation at these four points, f^ is discarded, and the 
function is evaluated at f^. At this point the right-hand extremal, e 2 » 
has been eliminated from the search which now inevitably proceeds to the 
left hand extremal at e.. . 

To find the second extremal, the function F is warped by writing 


HO = F(a) 


# 2H 

K = (ojj - j V 3 + a#i a>a * 

H 


(a6) 


), * 2N 

e = -(a - a L )[— f-~] + Ct*; a*>cx 




(at ) 
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where H is a positive integer, and a* is the location of the left 
hand extrenal. 

A typical relationship between £ and a is shown in figure (A 6) 
for the case K = 1. Differentiation of equation A7 with respect to a 
when IT * 1 results in 




; tt < a* 


(A8) 


Note that as a a* , 0 from both the left and right . At 

a = and at a * a^, £ =2* In the regions < a < a* and 

a* < a <cCjj £ varies parabolieally with a. Figure A7 illustrates 
these points. It can be seen that a region Act^ centered about a* 
transforms into a smaller region A£^ located in the neighborhood of 
£ = a*. On the other hand, a region Actg situated in the neighborhood 
of the upper search limit y naps into a wider region in. the neighborhood 

and a = ctjj are given 
by 2K; the greater the N, greater the warping becomes. 

The effeet of introducing, a moderate warping transformation on the 
function of figure is shown- in figure (A74* -It Can fee seen, from 
figure (At) that the region of influence of • . is reduced*,, and.* the 

region of influence of e^ is increased. On the warped surface search 
by sectioning commences with the evaluations of performance at to f^. 
Following these initial evaluations f^ is discarded (as opposed to the 


£ * a u . In general, the slopes at 




APPENDIX - B 


BUCLASP-2 Assumptions and Model 

This appendix is devoted to a discussion of some of the capabilities 
of BUCLASP-2 (A Computer Program for the Instability Analysis of Biaxially 
Loaded Composite Panels) as it pertains to the buckling analysis of the 
composite panels considered in the present work. This computer program 
(reference 10) is operational on the CDC-6600 computer. It is quite 
reliable and gives very good results for the buckling problems of compo- 
site panels. Some of the basic assumptions made in the analysis of 
HJCLASP-2 are as follows: - 

1. The panel members are orthotropic 

2. The material is linearly elastic 

3. Thin plate theory is employed 

U. Effects of prebuckling deformations are ignored 

5. Eccentricity effects are accounted for 

6. Exterior edges in planes normal to the prismatic direction 
are assumed to be simply supported. 

Support conditions at other boundaries are arbitrary. With the 
above assumptions an "exact” analysis of the whole panel is made. 

This analysis results in the prediction of Euler buckling modes, local 
buckling modes, or coupled Euler and, local modes. 

The user of BUCLASP-2 has to define the mathematical model of the 
panel under consideration. This mathematical model consists of three 
substructures, namely the start substructure, end substructure, and the 
repeat substructure. Figure B1 shows the cross sections of the three 
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substructures for the panel studied in this investigation. 

The results after using AESOP define the cross-sectional dimensions 
of the panel. These dimensions are used to find the buckling load using 


BU CLASP-2. 
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Figure 1.- Search processes. 
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Figure 2.- Optimization cycle. 
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Figure 3.- Panel to be optimized. 
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Figure 5.- Equivalent width. 
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Figure 6.- Weight strength plot for all-composite and all-aluminunVpaneiK, 
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Figure Al.- Search based on Golden section 
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Figure A2.- Response surface with two troughs. 
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Figure A4.- Non-convex response surface. 
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Figure A5.- Function with two extremes. 
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Figure Bl.- Mathematical model for BUCLASP2 
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Optimized design variables 
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Optimized design variables 
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ized design variables for all aluminum panel (L ■ 30 inches 




















